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The  linear  MHD  stability  of  a  confined  plasma  is  generally  studied  by 
means  of  energy  principles.  Up  to  date,  these  energy  principles  have  never 
been  justified  rigorously,  and  the  existence  of  a  solution  to  the  linearized 
equations  is  also  tacitly  assumed.  In  this  report,  based  upon  a  variational 
approach  Iwe  shall  first  establish  the  existence  and  uniqueness  of  a  gener¬ 
alized  solution  to  the  linearized  Lundquist  equations  for  a  toroidal  plasma 
confined  in  a  conducting  shell.  In  a  subsequent  report,  the  so-called  modi¬ 
fied  energy  principle,  which  includes  linear  energy  principle  as  a  special 
case,  will  be  justified  rigorously  and  a  solid  foundation  is  then  laid  for  the 
application  of  these  energy  principles  to  the  linear  MHD  stability  of  a  con¬ 
fined  plasma. 
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STABILITY  THEORY  OF  A  CONFINED  TOROIDAL  PLASMA 
PART  I.  EXISTENCE  AND  UNIQUENESS 

Peter  Laurence  and  M.  C.  Shen* 

§1.  Introduction. 

It  is  now  well  recognized  that  controlled  thermonuclear  fusion  may  offer 
the  ultimate  solution  to  the  world  energy  crisis  with  inexhaustible  resources 
for  years  to  come.  One  of  the  most  promising  methods  at  present  to  achieve 
the  goal  of  a  self-sustaining  thermonuclear  reaction  is  the  magnetic  confine¬ 
ment  of  a  plasma  inside  a  toroidal  vacuum  chamber  shielded  by  a  conducting 
shell,  such  as  a  tokomak.  General  surveys  about  the  progress  of  controlled 
thermonuclear  research  may  be  found  in  the  recent  articles  by  Grad  (1977)  and 
Furth  (1980),  among  others.  Needless  to  say,  physical  and  mathematical  prob¬ 
lems  abound  in  the  study  of  a  magnetically  confined  plasma,  such  as  equilib¬ 
rium,  stability,  wave  motion,  diffusion  and  adiabiatic  compression,  just  to 
name  a  few.  However,  at  present  the  problem  of  plasma  stability  still  remains 
the  central  issue  in  fusion  research,  especially  based  upon  magnetohydro¬ 
dynamic  (MHD )  models. 

Traditionally  the  standard  procedure  to  study  MHD  stability  of  a  plasma 
with  simple  geometry  is  to  linearize  the  governing  equations  about  a  given 
equilibrium  state  and  analyze  the  linearized  equations  via  normal  modes. 
Evidently  this  method  is  not  very  useful  wl  plicated  geometries  are 

considered.  To  avoid  this  difficulty,  Ber jt.ei -  al.  (1957)  developed  a 
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variational  approach  to  linear  MHP  stability,  and  formulated  the  so-called 
linear  energy  principle.  Its  advantage  lies  in  the  fact  that  to  determine 
stability  of  a  plasma  it  is  necessary  only  to  discover  whether  there  is  a 
perturbation  which  decreases  the  potential  energy  from  its  equilibrium 
state.  Since  then  a  vast  amount  of  literature  based  upon  this  technique  has 
appeared.  In  obtaining  the  necessary  condition  for  stability,  their  arguments 
rely  upon  the  assumption  that  the  eigenfunctions  of  a  certain  linear  operator 
must  form  a  complete  orthonormal  basis.  This  weakness  was  then  removed  by 
Laval  et  al.  (1965).  In  their  derivation  they  introduced  a  modified  energy 
principle  for  the  so  called  c-stability  of  a  confined  plasma,  which  includes 
the  linear  energy  principle  as  a  special  case.  Intuitively  speaking,  a  plasma 
is  called  o-stable  if  it  does  not  grow  faster  than  exp(ot),  where  <J  = 
and  T  is  soma  characteristic  time  needed  for  fusion.  This  concept,  indeed, 
is  more  optmistic  than  the  usual  exponential  stability,  since  in  fusion 
research  the  main  concern  is  whether  a  plasma  ean  be  confined  long  enough  to 
achieve  a  self-sustaining  thermonuclear  reaction.  On  the  other  hand,  many 
equilibria  regarded  as  unstsble  by  the  criterion  of  exponential  stability,  are 
in  fact  o-stable  and  may  be  found  practically  feasible  for  fusion 
research.  A  discussion  of  the  energy  principles  may  be  found  in  Greene  and 
Johnson  (1968),  Grossman  (1968),  and  Spies  (1976).  The  application  of  the 
o-stability  concept  to  various  plasma  configurations  was  further  expounded 
by  Goedbloed  and  Sakanaka  (1973,  I,  II). 

So  far  the  modified  energy  principle,  although  well  adopted  in  the  study 
of  MHD  stability,  has  naver  been  put  on  a  rigorous  basis  since  the  question  of 
the  existence  and  uniqueness  of  a  solution  to  the  linearized  MHO  equations  was 
never  asked  in  the  previous  proofs.  Without  a  precise  answer  to  the  question 
one  cannot  apply  the  energy  principle  with  complete  certainty.  The  first 
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attempt  at  the  existence  and  uniqueness  question  was  made  by  Rodriguez 
(1974).  Several  difficulties  however  were  found  in  his  proof.  Among  these, 
the  connectivity  of  the  vacuum  region  is  not  considered;  the  function  space  in 
which  the  Galerkin  scheme  is  introduced  is  not  complete;  the  scheme  does  not 
converge  in  the  norm  given  (which  in  fact  is  only  a  pseudo  norm).  Moreover, 
some  surface  terms  appearing  in  the  variational  formulation  are  discarded  by 
questionable  arguments.  In  fact,  these  terms  cause  grave  concern  in  obtaining 
estimates.  We,  on  the  other  hand,  will  make  certian  assumptions  about  these 
surface  terms  which  characterize  particular  equilibrium  states.  Indeed,  we 
doubt  that  solutions  would  exist  for  the  linearized  equations  in  the  absence 
of  such  an  assumption.  The  purpose  of  this  paper  is  to  give  a  correct  and 
rigorous  proof  of  the  existence  and  uniqueness  of  a  solution  to  the  linearized 
equations  governing  a  toroidal  plasma,  under  a  set  of  assumptions  imposed  upon 
an  equilibrium  state.  In  addition,  the  growth  rates  of  the  L  norms  of  both 
plasma  displacement  and  velocity  are  obtained.  In  a  subsequent  report.  Part  II 
we  shall  justify  the  modified  energy  principle  rigorously  and  a  solid  founda¬ 
tion  is  then  laid  for  its  application  to  problems  of  plasma  stability. 

In  the  following  ,  we  briefly  describe  the  contents  in  each  section.  In 
Section  2,  we  formulate  our  problem  based  upon  Lundquist  equations  for  magneto¬ 
hydrodynamics,  and  deal  with  the  linearization  of  these  equations.  Since  we 
consider  the  displacement  of  the  confined  plasma  from  its  equilibrium  config¬ 
uration,  it  is  natural  to  express  the  equations  in  terms  of  Lagrangian  coordi¬ 
nates.  Here  we  use  an  approach  which  is  essentially  that  of  Berstein  et  al. 
(1958)  so  that  the  final  results  may  be  checked  with  theirs.  In  the  vacuum 
region,  the  magnetic  field  is  both  curl  and  divergence  free.  We  follow  Berstein 
et  al.  and  introduce  a  vector  potential  instead  of  a  scalar  potential.  The 
latter  was  used  by  Lust  and  Martensen  (1960)  to  study  energy  principles. 


They  showed  that  the  two  approaches  are  essentially  equivalent  provided  cer¬ 
tain  cuts  are  made  in  the  vacuum  region  and  fluxes  through  these  cuts  are 
taken  into  consideration.  In  Section  3  a  mixed  boundary  and  initial  value 
problem  (MIBVP)  is  formulated  for  a  toroidal  plasma  confined  in  a  perfectly 
conducting  shell.  To  ensure  the  uniqueness  of  the  vector  potential  some  flux 
conditions  must  be  prescribed  (Blank,  Friedrichs  and  Grad,  1957),  and  a  set  of 
conditions  on  the  equilibrium  state  is  also  precisely  imposed.  In  Section  4 
the  Hodge  decomposition  theorem  is  used  to  define  function  spaces  which  play  a 
fundamental  role  in  this  work.  We  then  establish  the  continuity,  coerciveness 
and  symmetry  of  a  bilinear  form  in  Section  5.  In  Section  6  the  MIBVP  is  re¬ 
cast  in  the  form  of  an  evolutionary  variational  problem  (EVP)  by  means  of  the 
coercive  bilinear  form  (Lions  and  Magenes  1972).  We  then  construct  an  auxil¬ 
iary  EVP  (AEVP)  with  zero  flux  conditions,  whose  solution  is  established  by 
defining  Galerkin  approximations  and  showing  they  converge  to  a  solution  of 
AEVP.  The  solution  of  EVP  is  constructed  by  adding  to  the  solution  of  AEVP  a 
unique  solution  to  an  elliptic  system,  with  inhomogeneous  flux  conditions. 

The  uniqueness  proof  then  follows  easily  from  the  coerciveness  of  the  bilinear 
form  and  that  of  the  elliptic  system.  As  a  consequence  of  the  existence 
proof,  the  growth  rates  for  the  plasma  displacement  and  kinetic  energy  are 
also  obtained.  In  Section  7,  we  make  some  remarks  regarding  the  solution  of 
our  problem  and  the  relationship  between  a  solution  to  EVP  and  a  classical 
solution  to  MIBVP  is  discussed. 
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§2.  Linearization  of  Lundquist  Equations. 

We  consider  the  following  configuration  for  a  confined  plasma:  ft^  is  a 
doubly  connected  bounded  three  dimensional  region  surrounded  by  a  bounded, 
triply  connected  region  Thus  ft^  and  ftv  are  topologically  equivalent 

to  two  coaxial  tori  (Figure  1).  The  plasma,  a  highly  ionized  gas,  is 


contained  in  ft^.  The  vacuum  region  is  surrounded  by  a  conducting  shell 


In  the  plasma  region  ft  ,  (Bernstein  et  al,  1957) 

P 

P  ~~  =-Vp+uJXB, 
dt 

3B  _  „  x 

3t  -  -  V  x  E, 


V  x  B  =  pJ, 


(2.1  ) 
(2.2) 
(2.3) 


V  •  B  =  0, 


(2.4) 


§|  +  V(pv)  =  0,  (2.5) 


It  (PP~Y)  =  °' 


E  +  v  x  b  =  0. 


(2.6) 

(2.7) 


Here, 

u-81 

and  all  quantities  are  expressed  in  terms  in  Euler ian  coordinates,  v  is  the 
fluid  velocity,  B  the  magnetic  field,  J  the  current  density,  p  the 
scalar  fluid  pressure,  P  the  mass  density  and  Y  the  adiabiatic  constant, 
which  may  be  taken  to  be  equal  to  5/3. 


In  the  vacuum  region  ft^, 

7  x  B  =  o, 

V  •  B  =  0, 


(2.9) 

(2.10) 


the  vacuum  region  is  assumed  free  of  currents  or  charges,  and  the  fluid 
boundary  is  assumed  to  nowhere  intersect  the  outer  conducting  shell. 
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Furthermore,  the  two  regions  are  coupled  by  a  set  of  boundary  conditions 


where  [  ]  designates  the  saltus  of  a  quantity  across  the  fluid-vacuum 

boundary,  and  the  subscript  t  designates  the  tangential  component. 

On  the  conducting  shell 


Bn  =  0. 

In  addition,  since  is  a  multiply  connected  region,  if  the  Pre- 

Maxwell  equations  in  the  vacuum  region  were  solved  independently,  it  would  be 
necessary  to  prescribe  certain  "periods"  of  B  and/or  E  in  this  region  in 
order  to  ensure  uniqueness  of  solutions.  Physically,  the  system  is  linked  to 
the  interior  through  the  metal  shell.  However,  because  of  the  coupling  be¬ 
tween  and  the  situation  is  complicated,  and  we  will  not  return  to 

the  question  of  prescribing  periods  until  the  linearization  is  completed.  We 


also  note  that  a  detailed  discussion  of  the  formulation  of  the  boundary  condi¬ 


tions  may  be  found  in  Blank,  Friedrichs  and  Grad  (1957).  Also,  appropriate 
units  for  various  quantities  are  chosen  so  that  u  =  1. 

We  emphasize  that  the  non-linear  problem  defined  above  is  a  "free  bound¬ 
ary  problem".  U^(t)  and  hence  r^(t)  a  priori  unknown  must  be  determined  as 
part  of  the  solution.  Therefore,  we  wish  to  examine  the  effect  of  small  per¬ 
turbations  about  an  equilibrium  state  characterized  by  v  *  0.  This  equilib¬ 


rium  state  is  governed  by  the  equations: 


(2.14) 


In  Q  , 
P 


VD  =  v  X  B  x  B  , 
*0  0  0 


V  •  bq  =  0. 


(2. 15) 


In  a  , 
v 


V  x  b0  =  JQ,  (2.16) 


V  •  BQ  =  0, 


(2.17) 


with  the  boundary  conditions: 

1 

2 


B 


-]  =  0 

in 

r  , 

(2.18) 

po 

*  n 

on 

r  , 

(2.19) 

0 

po 

0 

on 

r  . 

(2.20) 

At  equilibrium,  we  may  assume  p  =  A(S)P  where  A(S)  is  a  function  of 
the  entropy  S.  Above  and  hereafter  a  superscript  p  or  v  will  designate 
plasma  or  vacuum  variables  respectively.  In  order  to  linearize  our  problem, 
we  reformulate  the  problem  in  reference  to  the  Lagrangian  coordi-nates : 

(x,y,z,t)  e  ap  x  r  — ■*  (X(),y0,z0,t)  e  aQ  x  R, 
where  (G(x,y,z,t)  =  (x0,y0,z0,t)  +  e5(xQ ,yQ ,zQt) . 


That  is,  we  identify  each  fluid  particle  by  its  position  in  the  equilibrium 
state  and  let  EC(x(J,y(j,z0,t)  be  the  displacement  of  the  fluid  particle 
xQ,y0,z0  from  its  equilibrium  state  xg*yo'zo  after  a  time  t.  He  first 
express  the  operator  V  in  Lagrangian  coordinates.  If  f(r,t)  is  given, 
where  r,t  are  the  Eulerian  variables  and 
f(r,t)  -  f(rQ,t), 

where  r  *  (x,y,z)  and  rQ  =  (xg,y0,z0),  then  it  follows  from  the  chain  rule 
that 

Vf (r,t)  -  V0f  •  Vr0. 
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Here  Vf  and  V  f  are  1*3  column  vectors,  and 


3  r  3r  3r 

7r  -  f  _ 2.  _ 2  — 2)  _ 

vrn  1  a„  »  a.,  '  a,  > 


3  x  3  x 
_ 0  _ 0 

3x  3y 


3y 

_0 

3x 

3zo 

L^T 


3y 

3z„ 


la* 

dz 

3z 

3z 


0 


Let 


r0  =  r  -  e^r0't}' 


and  express  C(r,t)  by 


C(r,t)  =  HrQ,t), 


S (rQ(r,t) ,t)  =  £(r,t). 


Then 


VrQ  =  X  -  eVQC  •  %, 

where  *  ^rQ  denote  the  usual  matrix  multiplication,  and  it  follows  that 


(X  +  eVQC)Vro  =  i. 

Formally,  we  may  assume  (I  +  eV^?)  1  exists  to  obtain 

(I  +  eV)_1  =  7r0' 


(2.21  ) 


(2.22) 


and  express 


(I  +  eVQS )_1  =  i  -  eVQ5  +  e2(VQC)2  +  •••  (-1  )n+1En(V^)n  + 


Without  attempting  to  justify  the  convergence  of  this  formula  or  its 
asymptotic  validity,  which  would  lead  us  to  consider  the  validity  of  the 
linearization,  we  henceforth  take 

VrQ  =  (I  +  tVo5)_1  =  I  -  e?o?,  (2.23) 

as  an  approximation. 
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We  gather  together  the  expressions  for  the  other  operations  involving  7 


as  follows.  To  first  order  in  E, 

V  •  v  =  Trace  (Vv) 

=  Trace  (7Qv  •  (I  +  E7^)”1) 

=  Trace  (VQv(I  -  e7QC) ) 

=  •  v  -  £  Trace(VQv  •  Vq£) 


and 


7  *  v  =  [ (I  +  •  Vq]  X  V 

=  [(I  -  *  7Q3  *  V 

=  [(VQ  -  £(VoS)t  •  vQn  x  v, 

where  t  denotes  the  transpose  of  the  matrix.  Furthermore,  note  that 

3e  *e  de  3* 

_+  u  .  V  =  —  =  ^  , 


where  superscripts  e  and  A  indicate  Eulerian  and  Lagrangian  coordinates 
respectively.  Thus  in  Lagrangian  coordiantes  (2.1)  to  (2.7)  become 

u  *i  . 

Pj^T  =  -VQp  •  [x  +  eVQ5)  +  ([(X  +  e%5)  ]  *  B)  *  B,  (2.24) 

|^  =  [(I  +  eV05)-1  7q]  x  (u  *  Bq )  +  (u  •  Vq)b  •  (I  +  eVoC)_1,  (2.25) 

|*T  =  -Yp  Trace  VqU  *  (I  +  e70^  ' *  (2.26) 


Here  the  second  and  third  equations  were  obtained  by  adding  u  •  Vb  and 
u  •  7p  respectively  to  both  sides  of  the  equation  (2.3),  and  (2.26)  is  the 
equation  which  results  from  combining  (2.6)  and  (2.7) 


3p 

3t  =  "u 


Vp  -  Yp7  •  u. 


Also  (2.5)  becomes,  in  Lagrangian  coordinates, 
3p 


=  P  Trace [V  •  u(I  +  e7QC)  1). 


(2.27) 
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We  now  let 


P  =  PQ  +  cp^t  p  =  pQ  +  cp^,  B  -  BQ  +  eB1#  u  ■>  Eu^ 
Since  zeroth  order  quantities  are  equilibrium  quantities,  we  obtain  the 
following  system 


'S 

— - - [-^(p„+Ep  )  (i-eV  £}]  +  — - - {i-eV  })fc  • 

po+epi  00  1  0  po+€pi  ° 

*  7o}  x  (Bo+eV  x  (VeV 

eB1 

{i-eV^^  •  Vq  x  {u1x(B0+EB1)}  + 

U  •  7o(B0+EBl)(I-E7o5) 

e*1 

= 

^(pjj+ep^  Trace {eV^  •  u^U-eV^)} 

Cp1 

(p0+e<V  Tracefe70»Ul  U-eVQZl] 

— 

_  CU1 

Keeping  only  first  order  terms  in  e,  that  is  equating  coefficients  of 

e  on  both  sides,  and  using  the  equilibrium  relation 

Vp  =  7  x  b  x  B 
Om)  0  0  0 

we  obtain 


— 

U1  ' 

I.  <v0p1-Op0»0„l, 

*  ^'VVV 

*  ^'VW 

p1 

-  Vo7o  •  »1 

3 

3t 

B1 

- 

?Q  X  (u1xB1 )  +  u  • 

Vo 

P1 

-  po7o  *  U1 

£  - 

-  °1 

.  (2.28) 


We  next  eliminate  all  quantities  not  directly  expressed  in  terms  of  u 
in  order  to  obtain  one  equation  for  u,  or  £,  which  can  be  solved  indepen¬ 
dently  and  then  used  to  compute  p.|f  Bj,  p^.  Alternately,  we  may  differentiate 
the  first  equation  in  (2.1)  with  respect  to  time,  assuming  that  spatial  and 
time  derivatives  commute  and  in  either  case  we  obtain  the  same  equation  but  in 

the  latter  case  for  u  rather  than  5.  The  equation  for  £  is 

°o  “a  5  ■  "Vo  •  {l  *  %  •  VW'-V  '  Bo 

ot 

+  Vo  X  (C  .  VoB0)  *  B0  +  Vo  X  B0  X  Vo  X  (CXB0) 

+  VQ  X  Bq  X  (5  •  VqBq)  -  (Vq5T  •  VQ)  X  bq  X  B0, 
which  reduces  to  the  equation  derived  for  P(5)  by  Bernstein  et  al  (1958) 

a2 

p0  — 2  *  ■  Voap070  *  5  +  C  *  Vo) 

3 1 

+  Vo  X  70  X  U;xb0)  *  B0  +  Vo  X  B0  x  70  X  (5^), 

provided  that 
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v  ( 5  •  V  p  ) 

0  0  0*0 


7  p 
0F0 


V  5  +  7  x  (5  •  7  b  )  *  B 
0  0  0  0  0 


+  vo  x  Bo  x  a  *  W  '  {V  *  vo)  x  Bo  x  Bo' 


the  derivation  of  which  is  given  in  the  appendix.  To  obtain  the  equation  for 
5  by  the  first  method,  we  integrate  from  0  to  t  the  last  four  equations 
in  our  system  to  obtain 

Pt(t)  =  -YpQ70  *  UU  +  Pl(0)  +  Ypo?o  •  5(0), 

B1(t)  =  V0  X  l^(t>  x  V  +  ^(t)  *  V0B0 

-  VQ  X  (5(0)  X  Bq)  -  5(0)  •  7qbq, 

Pl(t)  =  -P0VQ  .  5(t)  +  p^O)  -  pQV0  .  5(0). 

Assuming  that  our  system  passes  through  the  equilibrium  state  where  5=0, 

Pl  =  B1  =  0  at  some  time  tQ,  we  may  write 

Pl(t)  =  ~YPCV0  *  €(t), 

Bi(t)  =  Vo  x  (5(t)  x  Bo)  +  ^(t)  ’  VoBo' 

Pl<t)  =  -P0VQ  •  5 ( t ) , 


which,  when  we  plug  into  the  first  equation  of  (2.28),  together  with  the  last 
equation  in  (2.28),  again  leads  to  (2.29). 

The  linearization  in  the  vacuum  region  is  done  by  introducing  a  vector 
potential  A  there  such  that 

E  =  EQ+  EA,  B  =  BQ  +  eVQ  x  A, 

so  that  the  Pre-Maxwell  equations  there  are  taken  in  the  form 

VQ  x  V  A  =  0,  7q  *  A  =  0.  '  (2.29) 


If  we  are  dealing  with  the  equation  (2.29)  written  for  u  instead  of  5,  we 
retain  (2.29).  If  not,  we  again  assume  the  system  passes  through  the  equilib¬ 
rium  state  at  some  time  t  and  integrate  to  obtain 


7  *  V  *  a  =  0, 

0  0 


A  =  0. 
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For  the  boundary  conditions: 


(E  +  v  *  Bq  )  =0  or  n  x  (E  +  v  x  B)  =  0, 

become  to  first  order 

nQ  x  (A  +  C  x  Bq)  =  0  i.e.  n  x  a  =  (-hQ  •  £}B^, 

v  V 

or  in  integrated  form  n  x  a  =  -(n  •  £)b  ,  where  we  used  BQ  • 


furthermore 


2  2 

+  BP  BV  - 

P  +  —  =  ~  °n  Fp' 


may  be  written  as 


B  B 

(P0  +  +  <-J“  +  e2^2 


(_£  +  e ( 7  x  a))2, 


(BP)2 

(p0  -  £YP0?  *  C)  +  [— —  +  6(7  X  (£  X  B0)  +  £  •  VBq)] 

(B0)2  2 
=  (— ~  +  e(V  x  A) )  , 

which  if  we  retain  only  terms  of  the  first  order  in  t  reduces  to: 
-YPqV  •  C  +  BP  •  (5  *  Vb^  +  7  x  x  BP)) 


On  rv,  we  have  =  B^  •  (7  x  a  +  5  •  7B^) . 

Et=0,  or  nXE=0. 

Hereafter  we  will  omit  the  subscript  0  on  equilibrium  quantities. 
It  follows  that 

nXA=0  or  nXA=0  on  T  . 

v 

In  summary,  the  equations  derived  apply  either  to  pairs  (u^A)  without 
restriction  or  to  pairs  (C,A)  when  the  system  passes  throught  the 


equibl ibr ium  state. 
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§3.  Variational  Formulation  of  a  Mixed  Initial  Boundary  Value  Problem  (MIBVP). 


We  recapitulate  the  equations  in  §2  to  define  an  MIBVP. 

In  ft  , 

P 

p£  =  F(5)  =  V(YPV  •  5  +  5  •  7p)  +  !  *  B  *  V  *  15  »  B) 

+  V  x  V  x  (5  x  b)  x  B, 


in  ft  , 
v 

V  x  V  x  A  =  o,  (3.2) 

V  •  A  =  0,  (3.3) 

subject  to  the  following  conditions: 

(1)  Boundary  conditions. 

On  r  , 

P 


L((5,A))  =  0  =  -YpV  •  5  +  BP  •  <7  x  (^  x  BP)  +  5  •  Vb1*)  (3.4) 

-  BV  •  (7  x  a  +  C  •  7bV), 

n  x  a  =  -(n  •  ?)BV.  (3.5) 

On  r  , 
v 

n  x  A  =  0.  (3.6) 

(2)  Flux  condition. 

/  A(t)  •  n  ds  *  o(t),  o(t)  prescribed.  (3.7) 

r 

v 

(3)  Initial  conditions. 

5(0)  =  ||(0)  =  tQ,  MO)  -  Aq.  (3.8) 


In  a  natural  way,  we  shall  reformulate  the  equations  above  in  a  variational 
form,  which  is  more  amenable  to  analysis.  Its  equivalence  to  the  original 
formulation  will  be  discussed  later.  We  also  note  that  in  (3.8)  the  prescrip¬ 
tion  of  A( 0 )  is  not  necessary  since  it  can  be  determined  by  5(0). 
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We  shall  define  a  bilinear  form  on  the  product  space  a  *  ft  ,  denoted 

P  v 

by  a( (5, A) , (5, A) ) .  To  motivate  the  choice  of  this  bilinear  form  and  to  make 


clear  its  connection  with  the  original  problem,  we  proceed  as  follows: 

■v  -s. 

Given  a  pair  (C,A)  defined  on  £1^  *  £1^,  we  take  the  inner  product  of 
PC  =  F(C  ) , 


with  C,  integrate  over  £1  and  then  integrate  by  parts  to  obtain 

<PC,C>  ,  =  <F(C),C> 

t  (8  )  l/<ft  ) 

P  P 

»  /  {V(YpV  •  C+C  •  Vp)  +  V  *  B  *  ( V  x  C  x  B) 

£) 

P 

♦  ?  *  v  x  (?  »  bi  *  b}  •  {  dv 


-  /  “(YpV  •  {V  •  J  +  V  *  (5  *  B)  •  V  »  (?  *  B) 

£2 

P 

-£*Vxb*(Vx(£xB)-C*  V(5  •  Vp) 


+  /  (C  •  n){YpV  •  c  -  B  •  V  X  (t  X  b}, 

r 

p 

where  use  has  been  made  of  some  well-known  vector  identities.  In  a  similar 

>» 

fashion,  we  take  the  inner  product  in  a^  of  A  with  the  equation 

V  x  V  x  a  =  0, 

and  integrating  by  parts,  we  obtain 


0-+/  VxA*VxAdv+/  n  *  A  •  V  x  A  d v+J  n  *  A  •  V  x  A  dv. 

a  r  r 

V  p  V 

If  we  assume  that  the  pair  (C,A)  satisfies  the  boundary  conditions 

(3.5)  and  (3.6): 

n  x  a  =  -(n  •  C)BV  on  T  , 

P 

n  x  a  -  0  on  T  , 

v 

then  we  may  write 


0 


/  V  x  A 

a 


V  X  A  d  v  +  / 

r 

p 


(n 


C)BV  •  V  x  A  dv. 
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If  we  now  simply  add  this  result  to  the  right-hand  side  of  (3.7),  we  obtain 


<pl,l>  =  /  -tYp(V  •  5)(7  •  1) 

L2(fl  )  a 

p  p 

-  5  *  ?|5  •  Vp)  +  V  x  ((  x  b)  •  V  k  [5  *  B) 


-I*  V  x  b  x  ( 7  x  (5  x  B) ) }  dV  -  /  VxA'VxAdv 

Q 

V 

+  /  [?  •  n)[Yp7  *  £  -  BP  •  (7  X  (C  x  B5))  +  7  x  a  •  BV]ds. 

r 

P 

Let 

b( (£,A),(£fA) )  =  /  {Yp(7  •  £)(7  •  5)  -  5  •  7(5  •  7p) 

ft 

P 

+  7  x  (5  x  b)  *  (7  x  (5xB))-5*7xbx7x  (5  x  (5  x  B)  }dv 


+/  7  x  a  •  7  x  a  dv. 
n 

v 

and  it  follows  that 

<p5,l)  +  b( (5, A) , (5, A) )  = 

L  («  ) 

P 

/  (I  •  n)(Yp7  •  £  -  BP  •  (7  x  (5  x  bP)  +  7  x  a  •  BV)ds. 

r 

P 

Note  that  the  integrand  in  this  boundary  integral  differs  from  the  boundary 


condition  of  the  MIBVP  by 


So  by  defining 

2  2 
v  p 

a(  (5, A) ,  (5, A)  )  =  b((£,A),(£,A))  +  /  (£  •  n)£  •  7(2—  -  ~)ds,  (3.9) 

r  t  * 

p 

the  above  equation  becomes 
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<PC,C>  ,  +  a( (C,A) ,C,A) ) 

L  (Q  ) 

P 

=  /  (C  •  n){YpV  •  5  -  BP  *  x  x  bP)  +  V  x  a 

r 

p  v2  2 

+  €  *  V(3l-  -  sL>}ds 


BV] 


/  (C  •  n)L(  C  ,A )  ds . 

r 


(3.10) 


Thus  we  are  led  to  define  a  solution  of  the  evolutionary  variational  problem 
(EVP)  associated  with  the  MIBVP  as  a  pair  (C,A)  that  satisfies 


^  <pC,C>  , 

dt*  L  (fl  ) 

p 


+  a( (5, A), (5, A) )  =  0  , 


(3.11) 


for  all  (C/A)  in  a  space  W  to  be  defined  in  the  next  section.  After 
existence  and  uniqueness  of  a  solution  to  (3.11)  has  been  established,  we 
shall  return  to  the  question  of  the  equivalence  of  the  EVP  to  the  original 
problem. 

Remark : 

<v 

a( (C,A) , (C,A) las  defined  by  v3.9)  will  not  be  quite  right  for  our 

purposes,  however.  It  involves  a  term  of  the  form  /  C  '  7(£  •  Vp)<iv  which 

fl 

P 

will  not  necessarily  make  sense  for  (C,A)  in  W,  since  C  may  not  be  in 

H1 (fl  ).  Thus  the  form  of  a( (C,A) , (C,A) )  used  hereafter  is  obtained  by 

integrating  this  troublesome  term  by  parts  to  get,  from  (3.9), 

a((C,A),(C,A))  =  /  {Yp(V  •  0(7  •  C)  +  (7  •  C)(C  •  7p) 
fl 

P 


+  V  x  x  B)  •  V  x  (5  x  B) 


VxjjxVx(£xB)}dv 


2  2 
BV  BP 


/  7  x  A  •  7  x  a  dv  +  /  (C  •  n)C  •  7(~ - -  -  p)ds.  (3. 


12) 
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For  this  form  of  a( (4, A) , (C,A) ) ,  all  terras  as  seen  later  will  make 

sense  when  (£,A)  and  (5, A)  belong  to  W,  provided  that  the  measure 
2  2 
BV  BP 

dy  =  n  •  7(— - - — - p)  is  either  identically  zero  (no  surface  current)  or 

nonnegative. 

For  later  use,  we  have  to  make  precise  the  smoothness  conditions  on  the 
equilibrium  state.  Although  a  general  existence  theorem  (or  non-existence 
theorem)  for  a  three  dimensional  equilibrium  in  the  absence  of  symmetry  is 
still  lacking  and  there  is  good  reason  not  to  expect  one  except  under  special 
conditions,  the  assumptions  below  are  not  so  restrictive  in  comparison  with 
existing  exact  solutions.  We  assume  that 

p  e  h3(B  ),  Be  H3<n  ),  Be  H3(ft  ),  (3.13) 

P  p  V 

(for  the  definition  of  Hn(ft)  see  §4)  and  F  and  T  are  of  class  C4 

P  v 

and  have  the  cone  property.  It  then  follows  from  the  Sobolev  embedding 
theorems  (Lions  and  Magenes,  1972)  that 

p  e  ),  Be  c^n  ),  Be  c’fa  ).  (3.14) 

p  p  V 

From  the  trace  theorems  we  have  (Lions  and  Magenes,  1972) 

vP  e  h3/2(T  )  and  b  e  «5/2<r  )  n  H5/2(F  ), 

P  P  V 

which  imply  respectively 

p  e  c^r  >,  Be  c2<r  >  n  c2(F  ).  (3.15) 

P  P  V 

We  also  assume  that  n(x),  the  unit  normal  to  T  ,  can  be  extended  in 

P 

3  5/2 

U  as  a  function  in  H  (5)  ).  This  will  be  the  case  when  n(x)  e  H  (T  ) 

P  P  P 

because  of  the  inverse  trace  theorems.  Consequently  the  embedding  theorem 

implies  that  the  extension  of  n(x)  has  the  properties 

n(x) ,  V(n(x))  e  C°(fl  ).  (3.16) 

P 
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Furthermore,  the  equilibrium  is  assumed  to  describe  a  truly  sharp  boundary  in 


the  sense  that  there  exist  constants  c 
0  <  <  p ( x )  <  c2  <  “  in 

0  <  c,  <  P(x)  <  c  <  •  on 

3  4 

0  <  c,  <  |BV|  <  c,.  <  •  on 

D  D 


*  c4  '  c5 '  c6 


such  that: 

(3.17) 

(3.18) 

(3.19) 


Furthermore  we  assume  that  the  equilibrium  does  not  have  any  surface 

current,  or  that  when  there  is  a  surface  current  present, 

2  2 
P  v 

n  •  V(p  +  S- - 2__)  <  o  a.e.  on  T  .  (3.20) 

2  2  p 


He  conjecture  that  indeed  this  additional  assumption  is  necessary  in 
order  that  the  dynamical  equations  be  well  posed.  In  what  follows,  we  shall 
let 


(p  + 


(Bp)2 


(bV, 
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§4.  Function  Spaces 


The  object  of  this  section  is  to  introduce  various  function  spaces , 

culminating  in  the  introduction  of  a  Pre-Hilbert  space  W~  and  an  inner 

product  <•,*>  in  this  space.  We  shall  show  that  the  closure  of  the 

w~ 

space  W~  in  the  norm  II  ,  denoted  by  W,  preserves  certain  important 

w” 


functional  properties  to  justify  its  use  in  later  sections 
space  in  which  a  solution  to  the  EVP  (3.11)  is  sought. 

We  shall  denote  by  H^lS)  )  and  H^  (S)^)  respectively 
vector  valued  functions  f  =  (F..,f,,f,)  from  Q  and  SI 

I  £  J  p  V 


as  the  fundamental 

the  spaces  of 

3 

into  R  such 


that 

/  f2  d v  +  /  ( V f  .  ) 2 dv  <  -  Vi,  for  ), 

SI  1  SI  1  P 

P  P 

/  f2  dv  +  /  (Vf  J2dv  <  *  Vi,  for  H^S)  ). 

SI  1  SI  1 

v  v 

It  is  well  known  (Lions  and  Magenes,  1972)  that  the  space  H^SJ)  is  complete 
in  the  norm: 


■  f  ■  =  (I  If  J2  +  IVf  J2  }  ^  . 

H  (SI)  i  1  L  (SI)  1  LZ(Sl) 

Here  all  derivatives  are  taken  in  the  sense  of  distributions,  and  S)  denotes 


either  SI  or  SI  .  Similarly, 
p  v 


Hm(Sl)  =  {f  e  L2(Sl)|D0tfi  e  L2(S1)  for  |a|  <  m}  Vi  =  1,2,3, 


with  the  norm: 


m  =  (V  l  «Daul2  )  1/2  • 

H  (0)  i=1  |a|<m  L  (SI) 

Next,  to  define  a  space  we  first  recall  the  Hodge  decomposition 

theorem  (Friedrichs,  1955).  If 
v  e  L2  (SI), 

then  there  exists  v^  e  ,  V e  V^,  v^  e  such  that 
v  =  v,  +  v2  +  v3 
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or 


L  (ft)  =  v1  ®  v2  a  v3  , 

and  the  decomposition  is  orthogonal  and  thus  unique,  where 
v1  =  {V$  |4>  e  Hp(fl)}  , 

|  e  hJ  means  $  e  H1  and  <|>  =  0  on  Sft  the  boundary  of  ft, 
V2  =  {7  X  v|v  e  L2( rot )(«)}, 

L2(rot)  =  {v  e  l2(£1)|V  X  v  e  l2  (ft)}. 


(4.1  ) 


(4.2) 


(4.3) 


V3  =  { h | V  x  h  =  0,  7  •  h  =  0  in  ft, 


n  x  h  =  0  on  6ft} . 


(4.4) 


Remark:  If  V  e  L  (ft)  has  the  property  that 

7  •  v  =  0, 

then 

v  =  7  x  +  h 

where  ^  x  v2  e  ^2  and  h  e  V^.  This  is  a  consequence  of 

0=7*  v  =  A$  =*>  $  =  0,  since  e  e  H^. 

We  will  also  need  the  space 

L2(ft  )  =  (5|  /  p5c2dx  <  •}. 
p  P  ft 

P 

Note  that  given  the  assumption  (3.16) 


0  <  c_  <  P(x)  <  c  , 

3  —  —  A 

2 

the  space  t,  (ft  )  is  both  algebraically  and  topologically  equivalent  to 

p 

2  2 
L  (ft  ).  For  simplicity  we  shall  denote  the  scalar  product  in  Lp(ft^)  by 

2  2 

<«,»>2  p  ,  in  L  (0^)  by  <*»*>2  p»  and  in  L  hy  <*'*>2  v*  The 

corresponding  norms  are  denoted  by  I  l„  ,11  and  I  I  respec- 

2,p  2,p  2,v 

2  2 

tively.  If  it  is  clear  from  the  context,  we  shall  use  L  (ft^)  and  L  (ftv) 

to  denote  Lj^space  of  vector-valued  functions  or  scalar  functions.  We  shall 

also  use  P^fft^)  to  den°te  the  projection  of  the  set  W  of  pairs  (5, A)  on 

L2 (ft  ). 

P 
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We  now  define  the  space  W  : 


and 


{( C  ,A)  |C  e  hVj  ),  A  e  L2(rot)(3  )  n  V  (3  ) 
p  v  2  v 

V  x  V  x  a  =  0,  n  x  a  *  (-n  •  ?)BV  on  T  , 

P 

n  x  a  =  o  on  T  , 
v 


(4.5) 


W 


closure  of  W  in  the  inner  product 

<(C,A) , (C,A)>  =/  VxA'VxAd v  +  /  YpV  •  £V  •  £  dv 
3  o 

V  p 

+  /  V  X  (5  X  B)  •  v  X  (5  X  b) db  +  /  (C  *  n)(S  •  n)(-n  •  Vp)dx 

3  r 

p  p 


+  /  •  £  dv.  (4.6) 

3 

P 


To  establish  the  functional  properties  of  the  elements  in  W,  we  need  a 
fundamental  inequality  which  we  shall  first  prove  in  the  following 
Lemma:  If  A  e  L  (rot((3)  n  V2(S),  where  3  may  be  multiply-connected  with 

boundary  r ,  then 

I  Al  <  c{ I V  x  aI  +  In  *  A»  *.  }.  (4.7) 

2  2  H_/2<r) 

Proof : 

2 

By  duality.  We  will  show  that  when  A  e  L  (rot)  n  V2(3),  for  any 


f  e  L2(3),  a  =  (A^Aj/A-j),  f  =»  (fj,f2,f3,  then, 

|J  A  •  f  dv|  <  (IV  x  Al  +  In  *  Si  y.  }lfl_  . 

H"  *  <D  2 

For  the  proof,  first  note  that  we  can  reduce  to  the  case  where 

V2<3)  since  A  is  orthogonal  to  the  projections  of  f  into 

Let  f  be  so  given,  consider  the  auxiliary  problem, 

7  X  V  X  g  =  f 

in  3  r  =  ur, 

V  •  g  =  0  1  1 

n  *  g  =  0  on  T,  /  g  •  n  ds  =  0. 


f  e  L2  n 
V,  and  Vj. 


-23- 


i 

j 

» 

} 


Systems  of  this  kind  have  been  considered  in  the  literature  without  proofs  on 
specifications  on  f  by  Blank,  Friedrichs,  Grad  (1957),  by  Kress  (1970)  for 
the  case  of  infinitely  differentiable  f,  by  Solonikoff  (1978)  in  simply 
connected  domains  under  general  assumptions  on  f.  For  completeness  and 
simplicity  we  give  our  own  proof  of  the  existence  for  the  special  cases 
considered.  Let 


H(fl)  { 


g  C  L  (rot)fl  n  V^(fl),  n  *  g  =  0 


<g  •  n,1>  n  _ 1/  =0 

H  /2»  H  '2  (T) 


on  r , 

}  • 


H  is  a  Hilbert  space  for  the  inner  product 

<g»>H  -  <v  *  x  h>2  . 

This  follows  from  the  Friedrichs  inequality  of  the  first  kind  in  Friedrichs, 
(1955),  or  Werner  (1968). 

We  now  may  consider  the  variational  problem  to  find  g  such  that: 

/  V  X  g,  V  X  g  dv  =  /  f  •  g  dv,  V  g  e  H. 

n  n 

The  Riesz  representation  theorem  ensures  us  of  a  solution  g  e  H. 

Furthermore 

V  x  V  x  g  =  f. 


To  prove  this,  we  need  to  show  that 

/g»7xVxfdv=/  f  •  h  dv,  Vhe  c"(ft). 

Si  Si  0 

The  validity  of  (4.9)  follows  from  (4.8)  for  all  h  e  CQ  n  H(ft)  by 

SO  .o 

integration  by  parts.  To  prove  it  for  all  h  e  CQ(ft),  let  h  =  h  +  h  + 

3  13 

h  be  the  decomposition  of  Hodge  in  which  one  always  has  n  x  h  =  n  *  h 

*  0.  Since  heCg(fl),  n*h=0,  therefore  n  x  =  0,  and  h2  e  H(ft). 

Thus  since  h1  and  h3  contribute  nothing  to  either  side  of  (4.8)  because  of 

V  x  h*  *  V  x  “  0,  f  i  h\  holds  for  g  =  h  and  also  we  obtain  (4.9) 

again  after  integrating  by  parts. 
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Remark 


n  •  7  x  g  =  o  on  T. 

Note  that  n  •  Vx  is  a  first  order  tangential  differential  operator  on 

n  *  g,  and  therefore  since  n  x  g  =  0,  the  result  follows. 

Next  we  show  that  the  estimate 

«g«  ,  <  cifi  , 

H  («) 

holds.  Here  we  appeal  to  the  Friedrichs  inequalities  of  the  first  and  second 
kind  (Friedrichs,  1955) 

a)  Since  g  e  and  n  x  g  =  0,  we  have 

»g'5.(fl)  *  c,v  x  g,2,iJ  ’ 

b)  Since  n*Vxg=0,  n  x  g  =  o,  V  *  V  x  g  =  o, 

15  *  ’V(!»  ‘  ',v  x  '  *  ^2,a- 

The  result  now  follows  since 

7  x  (-5—  g)  =  — ^-j-(7  x  g)  in  D' . 

3x  3x 

We  now  return  to  the  proof  of  the  Lemma, 

|/a*  f  dv|  <  |  /  A  *  V  x  V  x  g  dv| 

n  n 

<  |  /  V  X  A  •  V  »  9  dv|  +  |  J  n  *  A  •  v  g  dv| 

n  r 


<  HV  x  A  I  x  9I  +  In  X  Al  v  IV  x  gl  . 

2'Q  2'°  H  ^  (D  h  V2(D 

<  17  x  AS  q17  x  gl?  n  +  c,n  *  A*  -V  ,9'  2  * 

2,“  2,“  h  /2(D  h2(Q) 

where  for  the  last  inequality  we  have  used  the  trace  theorem,  lul 
cIu*h,^qj,  (Lion  and  Magenes,  1972)  given  the  smoothness  of  T. 
Thus  since. 


H  1/2  (T) 


»gi  ,  <  »f«2  8. 

HZ(fl)  Z' 


and  since  obviously 


"7  *  o  «  cl911  , 

H2(fl) 
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there  follows 


|  /  A  *  f  dv |  <  c{lV  x  A|  +  In  x  Al  ,  }«fl, 

0  2,Q  H_/2(r)  2'^ 


This  proves  the  Lemma. 
Remark . 


Note  that  if  (-n  •  V  p)  >  c  >  0  a.e.  on  T  ,  the  proof  would  have 

P 

been  considerably  simpler,  since  then  the  more  elementary  estimate, 


«Al2  Q  <  c{lV  x  A«2  +  In  x  Al 

'  L  (T) 


} 


may  be  used. 
Theorem  4.1. 


I. 

II. 

III. 

IV. 


is  a  Hilbert  space  of  pairs  (C,A)  satisfying 

C  e  L2(fl  ),  v  x  (C  X  B)  e  L2(fl  ),  v  •  c  e  L2(ft  ). 
p  p  p 

a  e  L2(ft  ),  v  x  a  e  L2(ft  ),  V  x  v  x  A  =  o. 

V  V 

n  x  a  =  (-n  •  C)BV  on  T,  nxA=Ot  on  T. 

P  v 

W  is  separable. 


Proof  j 

I .  We  show  C  e  L2  ( ft  )  =  L2  ( ft  ) . 

P  P  P 

For,  assume  that  (Cm,Am)  e  W  is  a  Cauchy  sequence  in  the  W~  norm, 

then  by  (4.6)  1C  I  „  is  Cauchy, 
m  2 ,  P 

aC  e  l2(0  ),  icm  -  CI  ♦  o, 

P  2  ,p 

2  2 
since  L  (0  )  is  complete.  It  is  also  easy  to  see  that  C  e  L  (ft  ). 

P  P 

Next  we  show  that  V  x  (£  x  b)  6  L2(ft  ). 

P 

If  (Cm,Am)  e  W  is  a  Cauchy  sequence,  then  by  (4.6) 

IV  x  am  x  B) I 

2»P 

is  Cauchy,  and 

3G  e  L2(ft  )  such  that  IV  x  (5m  x  B)  -  Gl  ♦  0. 

P  2 ,  p 


and  in 


Also  3£m  e  L2(ft  )  such  that  £m  x  B  ♦  £  x  B  in  I  l_  . 

P  P  2,P 

particular  D’.  Hence 

7  x  (Cm  x  B)  7  x  (5  x  b). 

So  by  the  uniqueness  theorem  for  D1, 


7  x  (C  x  B)  =  G  e  L  (ft  ). 

P 

2 

Finally  we  show  7  •  £  e  L  (ft^). 


We  make  use  of  the  assumption  (3.15)  on  the  pressure.  If  (£m,Am)  is  a 

ID  2 

Cauchy  sequence  in  I  I  ,  then  Y  7  •  £  is  Cauchy  in  L  (ft  ).  Hence 

W  p  P 


is  Cauchy  since  p  is  bounded  below.  Thus,  by  the  same  reasoning  as 


above,  7  •  £,  taken  in  the  sense  of  distribution,  is  in  L  (ft  ). 

P 

II.  A  e  L2(ft  )  and  7  x  a  e  L2(ft  ),  7  x  7  x  A  =  o. 

v  p 

Assume  (£m,Am)  converges  in  I  1^,  then  in  particular  there  exists 

m  o  5 

W  e  L  (ft  )  such  that  7xA  ♦  W  in  L  (ft  )  and  there  exists  £  e  L  (ft  ) 
V  v  p 

such  that  £m  +  £  in  L2 (ft  ) ,  7  •  +  7  •  £  in  L2(ft  ).  These  last  two 

P  P 

results,  together  with  the  inequality  (Teraam,  1974): 


imply 


•n  •  ni  !  <  c{in«  +  17  •  ni  } 

h"  /2  ( r )  l  (ft  )  u  ( ft  ) 

p  P 

£™  •  n  ♦  £  •  n  in  H  ^2  ( T  ) . 


(4.10) 


Therefore  since  by  assumption  Bv  is  in  C^fT^), 

(£  •  n)BV  e  H~  (T  )  , 

2  -1/, 

for  ♦  e  L  (T  ),  dense  in  H  2  (T  ),  we  have 
P  P 


n  B  1/_  _  1/n 

h  /2  x  h  /2  ( r  )  r 


/  (?  *  si  BV)  •  ♦  ds 


/  (£  •  m)  BV  •  $  ds. 

r 


Note  that 


/  (5  •  m)  BV 

r 


$  ds  <  15  *  ml  |.  IB 

H_  '2  (T  ) 

P 


♦  ■  v, 

H  /2  ( T) 


<  cl£  •  ml 


-V,  '♦*  V, 

H  '2  (T  )  H  '2  (T  ) 
P  P 


-27- 


Hence  elso 


(Cm*  b)BV  ♦  (£  •  m)BV  in  H  2  (T  ). 

P 

Now  since, 

n  x  Am  »  (£m  *  n)BV, 

n  x  A™  is  Cauchy  in  H  ^2  ( T  ) . 

P 

Therefore  the  statement  A™  +  A  in  L2(£2  )  follows  by  the  fundamental 

v 


inequality  (4.7) 

I Al _  <  c{ 17  x  Al  +  In  x  aI 


2,v 


2,  v 


h“  1/2  ( r  > 


V  A  e  L  (rot(fl  )  n  V,(8  ). 

v  2  v 


Note  that 


In  x  aI  i,  =  In  *  aI  y.  , 

H“  /2  ( T )  h“  /2  (T  ) 

P 

since  In  x  aI  i,  ■  0  because  of  boundary  condition  (2.6). 

h“  /2  <r  j 

V 

III.  n  x  a  =»  (-n  •  C)BV  on  r  ,  n  x  a  =  0  on  T. 

P  v 

a)  (-n  •  Vp)  >  c  >  0  a.e.  on  T,  (n*C)m  is  a  Cauchy  sequence  in 

P 

l? (T  )  directly  from  our  definition  of  the  norm  l(€,A)l  ,  thus 

P  m  2  w_ 

(n  x  a)  is  a  Cauchy  sequence  in  L  (T^)  since  Bv  is  smooth. 

Furthermore  (n  x  A)m  is  Cauchy  in  L2(T  )  since  n  x  a"1  =  0  on  T  . 

v  v 

Passing  to  the  limit  in  the  boundary  condition,  we  have 

n  x  a  *  (-n  •  6)BV  in  L2(T  ) 

P 

nXA®0  in  L2(T). 

v 

b)  When  (-n  •  Vp)  =  0  on  T  on  a  set  of  non  zero  measure,  we  make  use 

P 

of  the  estimates  (4.7)  and  (4.10)  to  conclude 

(-n  •  C)BV  -  n  x  a  in  H~^2  (T  ) 

P 

n  x  a  -  0  in  H~  ^  ( T  ) . 

V 
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VI.  Finally,  we  show  that  W  is  separable. 

This  is  so  because  W  can  be  embedded  in  a  closed  subspace  of  a  product 
2 

of  L  spaces  in  the  following  way: 

W  **•  L 2  (2  >  X  L2(£2  )  X  L2(fl  )  x  L2(fl  ) 

p  p  p  V 

(5, A)  C — *■  (£,YpV  •  £,V  x  (5  x  B),V  x  a).  (4.11) 

Here,  we  take  as  a  natural  norm  on  the  product  space,  the  one,  whose  square  is 

2  2 

the  sum  of  the  squares  of  the  L  norms  on  each  L  space. 

Thus  by  virtue  of  the  definition  of  •  I  (4.6),  this  mapping  (4.11)  is 
clearly  an  isometry  of  W  onto  its  range,  hence  the  range  is  closed.  Thus 
W  is  separable  since  any  closed  subspace  of  a  Hilbert  space  is  separable. 

Another  important  fact  necessary  to  understand  the  structure  of  the 
space  W  is  the  following. 

The  vacuum  member  A  of  a  given  element  (£,A)  e  W  is  uniquely  deter¬ 
mined  by  the  boundary  values  £  •  n  on  T  .  He  formulate  this  statement  as 

P 

Theorem  4.2. 

There  exists  a  unique  solution  of  the  system: 

V  x  V  x  a  =  0,  in  SI, 

v 

V  •  A  =  0,  in  S2 

v  ' 

subject  to 

n  x  a  =  -(n  •  £ )BV  on  T  ,  (4.12) 

P 

n  x  a  =  0  on  r  , 

v 

and 

/  A  •  n  ds  =  0. 

r 

v 
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Also,  this  solution  A  is  in  V  •  In  particular,  note  that  the  only  solution 
of 

V  x  V  x  a  =  0 ,  V  •  A  =  0,  (4.13) 

n  x  a  =  0  on  T  u  T  ,  /  A  •  n  ds  =  0  , 

p  v  r 

V 


is  the  solution  A  =  0.  The  proof  of  this  proposition  is  given  in  Kress 
(1970). 


To  conclude  this  section,  we  define  three  more  spaces  for  later  use: 
H  :  closure  of  W  (or  W~)  in  the  norm: 


-+ 


1  ,  =  {l£l 

l  +  iA» 

H 

2  ,P 

W+  and 

H+.  Let 

A 

{(C,a)K 

e  Pn  (w  ) 

(4.14) 


where 

A  e  Pfl  (W-),  A*  e  v3 
v 

then 


w+  - 

closure 

of 

w“+ 

in 

the 

W  norm. 

H+  = 

closure 

of 

w“+ 

in 

the 

H  norm. 

(4.15) 


(4.16) 

(4.17) 
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§5.  Continuity,  coerciveness  and  symmetry  of  a((£,A),  (£,A)). 

The  expression  for  a( (C,A), (C,A) )  defined  in  §3  is  given  by 

a((£,A),(OA))  =/  5  »  A  •  V  *  J  dv 

ii 

V 

+  /  {V  x  (5  X  B)  •  v  X  (C  X  B)  +  (V  •  C)(C  •  Vp) 
ft 

P 

-  C  •  V  x  B  X  (V  X  (C  X  B ) )  +  Yp(7  •  0(7  •  Odv 


-  /  (€  •  n)<5  *  n)n  •  Vpds  ,  (5.1) 

r 

P 

where 


We  first  prove  that  a( ( C,A) , (£,A) )  is  a  continuous  bilinear  form  in 
W  x  w. 


Theorem  5.1. 

»v  <v 

For  any  (€,A),(Oa))  e  w,  there  exists  c  >  0,  a  constant,  such  that 
|a(  (C ,A) ,  (Oa)  >  |  <  cK5,A)lwKU)lw  •  (5.2) 

Remark : 


Hereafter  we  shall  use  c  as  a  generic  positive  constant. 

Proof:  The  integrals  on  the  right  hand  side  of  (5.1)  are  estimated  as 
follows: 

/  V  x  a  -  V  x  dv  <  (|  |V  *  A|2dv}{/  |7  x  A|2dv}  ^ 2  .  (5.3) 

n  n  n 

V  V  V 

/  YpV  •  57  •  X  dv  <  {/  Yp{V  •  C)2dv}  1/2  {/  Yp(?  •  02dv}  ^ 2  .  (5.4) 

ft  an 

p  p  p 

/  V  x  (5  x  B)  •  v  X  (C  x  B)dv 

a 

p 

<  If  |V  X  (5  x  B) |2dv}2{/  |7  X  (X  X  B) I2dv}  ^  .  (5.5) 

n  n 

p  p 
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where 


/  |{  •  V  x  b  *  |V  *  (5  x  B) | dv 

a 

P 

<  (/  »?2a.)  V*  (/  l*  «  <t  »  b>|2<m  v2 

a  n  p 

p  p 

2 

<  {/  PS  dv}  ^  {c  J  |  V  x  ( £  x  B)|2dv}  ^  . 

ft  ft 

P  P 

/  (V  •  C)(£  •  V  ) dv 

ft  P 

P 

<  (/  Yp(V  *  C)2dv)(/  Pt2dvj 

ft  ft  Tpp 

P  P 

<  c(/  Yp(V  •  C)2d v) ( /  pC2dv). 

ft  ft 

P  P 

/  (C  *  m)(C  •  n)(-n  •  Vp)ds 

r 

P 

<  /  (C  •  m)2(-n  •  Vp)ds  /  (C  *  m)2(-n  •  Vp)ds, 

r  r 

P  P 


&•  6* 


on  r  . 
p 


(5.6) 


(5.7) 


(5.8) 


(5.9) 


(5.10) 


-n  •  Vp  >  0 

In  (5.3)  to  (5.9)  the  essential  tool  is  the  Schwartz  inequality.  In 
(5.4)  we  apply  it  after  writing  Yp  *  ^Yp  /y p.  In  (5.6)  and  (5.8)  we  first 
multiply  and  divide  by  /p.  In  (5.8)  we  multiply  and  divide  by  /yp  .  In 


(5.10)  we  use  the  Schwartz  inequality  after  writing  -n  •  Vp  =  /-n  •  Vp 


/-n  •  Vp  . 

Thus  (5.3)  through  (5.10)  combine  to  yield 

|a((S,A),(?,A)|  <  c(l/Yp  V  •  Cl  l/yp  V  •  «l 

2  *p  2,p 


+  IV  (C  X  B) I  IV  x  (C  x  B) I  +  l/Yp  V  .  Cl  I ci _ 

x  2,p  2,p  *  2,p  2,p 

+  IV  x  (C  x  B) 1 2  Ul  +  {/  <C  •  n)2(-n  •  Vp)ds }{/  (C  •  m)2(-n  •  Vp)ds}. 

2#P  2,p  r  r 

P  P 
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Therefore  using  Schwarz's  inequality 


l  a.b.  <  a2)  1/2  (SI  b2)  , 

U  1  X  1  1 

we  obtain, 

|a((£,A),(£,A)|  <  c{«/y5  V  •  +  |7  x  (C  x  B)|2  ♦  lei2 

2,p  2,p  i, P 

+  /  (C  •  m)  2ds}  { l/7p  7  .  hi  +  17  X  (X  X  B)|2 
P  2»P  2 ,p 

P 

+  i'll2  +  /  (?  •  m)2ds}  <  cl  (€,A)  I  I  (?,A)  I  . 

2  ,  p  p  W  W 

p 

This  completes  the  proof  of  the  theorem. 

Next  we  show  that  under  the  assumptions  (3.13)  to  (3.18), 
a( (C,A) , (S,A) )  is  coercive. 

Theorem  5.2. 

There  exist  positive  constants  6  and  X  such  that: 

a((C,A),(C,A))  +  Xl£l2  >  «MS,A)I2  ,  (5.11) 

L2(S1  ) 

where  P 

a((£,A),(5,A))  =  /  (Yp(V  •  i)2  +  [7  x  (£  x  B)]2 

a 

P 

+  (7  •  £)(£  •  7p)  -5»7xBx7x(5x  B) }dv 


I  /  (7  x  A)2dv  -  /  (5  •  n)2n  •  7p  ds  . 

SJ 

v 

An  integration  by  parts  of  the  volume  terms,  first  shown  by  Kruskal  et 

al.  (1958),  justifies  the  identity  (Spies,  1974): 

/  (Yp(V  •  5)2  +  [7  x  (C  x  B) ] 2  +  (7  •  5)(£  •  VP) 

fl 

P 

-£*7xBx7x(£xB) 

-  /  [7  x  ((  x  B)  +  n  •  57  x  B  x  nj2  +  yp(7  •  £)2 

fl 

P  2 

-  2(V  x  B)  x  n  •  (B  .  7n)(£  •  n)  dv  . 
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Here  the  only  possibly  negative  quantity  is 


-  /  2(V  x  b)  x  n  •  (B  •  ^n)(£  *  n)2dv, 

ft 

P 

and 


Letting 

2(V  x  b)  x  n  •  (b  •  Vn)(x)  =  M(x), 


we  may  find  a  constant  e  such  that 

-  /  M(5  •  n)2dv  +  e  /  p£2dv  >  0  , 

ft  ft 

P  p 

V5,  such  that  (C,A)  e  W. 

Since 

C2  =  (n  •  £)2  +  (n  x  (n  x  C))2  ' 

P(€  •  n)2  <  PC2  , 


it  suffices  to  let 


»  B  «  n  ■  B  ■  Vnf 
P 


(5.12) 


(5.13) 


So  by  integration  by  parts,  we  obtain 

a( (C,A) , (£,A) )  +  e«CI2  p  >  0  ,  (5.14) 

V(?,A)  e  W. 

(5.14)  will  prove  a  useful  result  when  we  analyze  the  stability  of  solutions 

to  the  MIBVPj  however,  it  is  not  sufficient  to  establish  coerciveness  as  it  is 

2 

not  possible  in  this  approach  to  obtain  5l(£,A)lw  on  the  R.H.S.  of  (5.11). 
This  can,  however,  be  done  in  the  following  ways 
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I /  {  .  V  X  B  *  V  x  (5  x  B]dv| 

ft 

p 

<  /  |5| |V  x  B| |V  x  (£  x  B) |dv 

n 

p 

<  e  J  |5| |V  X  (5  X  B)|dv 

n 

p 

<  e,{/  K  |2>  1/2  {/  |V  x  (5  x  B)>2}  1/2 

1  n  a 

p  p 

<  e,  /  K|2dv  +  f.  /  |V  X  (5  X  B)|2dv 

1  4f1  ft  1  ft 

P  P 

<ei  {5i;  k,2,p  +  Vv  x  axB),2,P}' 


where  we  have  used  the  inequality 

ab  <  j;  a*4  +  f*>  , 

4f 

and 

e  *  IV  x  bI  „  ,  e2  =  [gib  P(x>]  1  . 

L 

In  the  same  way,  we  obtain  that,  for  arbitrary  f2  » 

e  I5I2  2 

|  /  (V  •  £)(5  •  VP)dv|  <  e  ( . 24T  2~^  +  f2e4l  ^  *  *  £'2/p}  , 

ft  2 

P 

where 

e  =  iVpl  „  ,  e  -  [glbYp]"1. 

J  L 

The  surface  contribution  in  the  expression  for  a( (5, A), (5, A))  is 

2  2 

-  /  (5  •  n)2n  •  V(p  +  5- - ®_)<jv 

r  2  2 

P 

which  is  nonnegative,  by  assumption. 
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We  now  piece  together  the  above  estimates  and  find 

e  e  e  e 

a((S,A).(£,A>)  +  +  ll'^l  22p 

>  (1  -  e  f  )IV  x  (£  x  B )  1 2  +  (1  -  e  e  f  ) I /ip  V 

11  4f  P  J  4  <c 

+  /  (S  •  n)2(-n  •  Vp)ds  +  «v  X  Al2  +  . 

p  4|P 

P 

By  choosing  and  f 2  for  instance  so  as  to  satisfy 


e1f1  =  V->  e3e4f2  =  » ( 5. 15  ) 

we  thus  obtain  from  (5.15) 

a(  (€,a)  ,  (€,a)  )  +  Al?l^  >  <51(5, A)l2 

4  f  P  W 


with 


2f, 


2f„ 


6  =  V2  , 

1  ""2 

Finally  we  establish  the  symmetry  of  a( (t,A) , (£,A) ) . 


Theorem  5.3. 

For  any  pair  (5 /A)  e  W+,  (£,A)  e  W+  , 

a( (C,A) ,  (S,A))  =  a((€,A),(£,A)). 


£1 


2 

2,  P 


(5.16) 


Proof  s 

From  (5.1),  for  (S,A)  e  W_+,  (C,A)  e  w”+ 
a((?,A),(?,A))  -  a( (C,A) , ( C,A) )  = 

/  U  •  V(C  •  Vp)  -  5  •  V(5  •  Vp)  +  £  •  V  x  b  x  V  x  (£  x  b) 

n 

p 

-£*VxbxVx{£x  B) }dv.  (5.17) 


If  we  can  show  that  this  expression  is  zero,  then  since  W-+  is  dense  in 
W+,  by  passing  to  the  limit  in  pairs  (5m,Am)  e  W  +  W+  (£,A)  e  W+, 

(5  ,A  )  e  W  +  (C,A)  e  W  ,  the  proof  of  symmetry  of  W  will  be 
complete . 

We  will  show  that  the  expression  above  can  ?  expressed  as  a  divergence 
which  in  turn  can  be  converted  into  a  surface  term  which  vanishes. 
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Now  we  establish 


/  (-5  *  7(£  •  v  )  +  £  •  7(£  •  Vp)  +5«7xBx7x(£xB) 

ft  p 

P 

-5*7xBx7x(£xB) }dv 

-  +  /  7  •  {B  •  (5  *  5 )  J  }  dv  . 

ft 

P 

For  this,  we  show  the  following  identity  holds: 

B(5,5)  =  -5  •  7(5  •  7P)  +  6  .  7(5  •  7P) 
+C»Vxbx7x(5xB)-5.7xBx7x(5xB) 

-  7  •  (b  •  (5  *  5)j)  =o  (5. 18) 

for  5,5  e  P^  (W  +  ),  the  set  of  all  first  components  of  pairs  (5, A)  e  w~+. 

P 

Note  that 

b<5  +  5\5  +  5')  =  B<5,5)  +  B(5',5)  +  B(5,5')  +  b(5',5').  (5.19) 

If  the  pair  (5,5),  <5\5),  (5,5'),  <5',5')  satisfy  the  identity,  then  the 

<N  »%» 

pairs  (5  +  5', 5  +  €')  also  do.  Also  note  that  if  either  the  first  or  second 

<v 

member  in  the  pair  (5,5)  is  zero,  the  identity  is  trivially  satisfied. 

To  establish  (5.18),  we  need  the  Lemmas  below. 

Lemma  5.1. 

B(5,5)  =  0  implies  B(a(x)5,5)  =  0  , 

V  scalar  a(x).  (5.20) 

Proof : 

(5  •  7a) (5  •  7P)  -  7  x  b  .  5  x  <<7a)  x  (5  x  b)) 

=  (7  x  B  •  7a) (B  •  5  x  5)  . 
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That  this  last  equation  is  indeed  an  identity  is  immediate  from  the  relations 
(£  *  Va)<£  •  7p) 

=  (5  •  7a)(£  •  7  x  b  x  B)  -  7  x  b  •  £  x  (7a  x  (£  x  b)) 

=  7  x  b  •  (£  •  (£  x  B)7a  -  (?  •  7a)(£  x  b)) 

=  (7  x  b  •  7a)(£  •  (5  x  b)  +  ( £  •  7a)  •  (£  •  7  x  b  x  b) 

=  -(7  x  b  •  7a) (b  •  (£  x  X)  +  (£  *  7a)(5  •  7  x  b  *  B). 

Lemma  5.2. 

The  identity  (5.18)  holds  for  the  Cartesian  basis  vectors 
£  =  e^  ,  £  =  e_.  ,  i,j  =  1,2,3. 

The  proof  of  Lemma  5.2  is  straightforward  but  tedious  and  we  omit  it  here. 
Lemmas  5.1  and  5.2  together  with  (5.19)  combined  with  lineariy  and  density 
of  vT+  in  W+  establish  Theorem  5.3  . 
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§6.  Existence  of  a  Solution  to  EVP. 

In  this  section,  we  establish  the  existence  of  a  solution  to  the 
following  EVP: 

In  the  interval  (0,T],  T  >  0,  find  (5,A)(t)  e  W  such  that 
2  ^  * 

~  <C(t) ,  £>_  +  a((S,A)(t M?,A)>  =  0  , 

dt2  2,P 

V(C,A)  e  W+  (6.1 ) 

subject  to  the  initial  conditions 

(€(0 ) ,A( 0 ) »  =  ( ,Aq )  e  W+  ,  (6.2) 

||(0)  =  lQ  e  Pfl  (H),  (6.3) 

p 

and  the  flux  condition 

/  A(t)  •  n  ds  =  0(t),  (6.4) 

r 

v 

where  o( t)  is  a  given  continuous  function  of  t  on  [0,T).  The  problem 
will  be  divided  into  two  stages.  Much  of  the  work  consists  in  defining  a 
solution  to  an  auxiliary  evolutionary  variational  problem  (AEVP)  and  proving 
the  existence  and  uniqueness  of  a  solution  to  the  AEVP.  Then  we  shall  add  to 
the  vacuum  part  of  the  solution  of  this  auxiliary  problem  an  appropriate  solu¬ 
tion  of  an  elliptic  system  of  equations  in  order  to  obtain  the  actual  solution 
to  the  EVP. 

§6.1.  An  Auxiliary  Evolutionary  Variationaj.  Problem  (AEVP). 

The  AEVP  is  formulated  as  follows: 

On  [0,T],  find  (£,A)(t)  e  C([0,«):»)  with  lit)  e  C([0,»):Pfl  (H)) 

P 

such  that 

2 

~  <5(t),?>  +  a((S,A)(t),(t,A))  =  0  , 

dt 

V(?,A)  e  W,  (6.5) 
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subject  to  the  initial  conditions 


(C(O).A(O))  =  (C0,A0  -  Aq )  e  w. 


(6.6) 


ff(0>  =  ^0  e  p«  (H+)  =  pfl  (H)' 

p  p 


(6.7) 


where  specification  of  A(0)  is  not  required,  ano  A^  is  an  element  of 
satisfying 

/  A'  •  n  ds  =  3(0),  (6.8) 

r 


and  the  flux  condition 

/  A  •  n  ds  =  0. 

r 


(6.9) 


Note  that  this  flux  condition  actually  follows  (cf.  §4)  from  A  e  (W),  but 

v 

we  state  it  explicitly  here  for  the  sake  of  convenience. 

To  show  the  existence  of  a  solution  to  the  AEVP,  we  shall  define  a 
sequence  of  Galerkin  approximations,  and  then  use  the  coerciveness  of 
a( (5, A) , (C,A) )  to  derive  an  energy  inequality.  This  energy  inequality  will 
then  be  used  to  show  that  the  solutions  are  strongly  bounded  in  a  suitable 
space,  and  hence  converge  weakly.  Finally  we  shall  show  that  the  limit  is 
indeed  a  solution  of  AEVP.  In  order  to  define  the  Galerkin  approximations, 
we  begin  by  proving  the  following 

Lemma  6.1:  There  exists  a  sequence  of  elements  8*  “  that  satisfy 

the  following  conditions: 

(1)  {0  constitutes  a  basis  in  H  and 

<0>?>2,P  -  *ij  ' 
i. 

(2)  {0  form  a  complete  set  in  W,  i.e.  finite  linear  combina¬ 

tions  of  elements  of  the  sequence  are  dense  in  W. 
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r 

» 


Proofs  Since  W  is  separable,  there  exists  a  sequence  of  linearly  indepen¬ 
dent  vectors  {a*}  which  form  a  basis  in  W.  We  now  use  the  Gram  Schmidt 
orthogonal izat ion  process  to  othonormalize  the  first  components  ct  in  the 

norm  I  S as  follows.  Let 

2,P 

B1  .  aVlaJl 

IT(U  ) 

P  p 

B2  =  {a2  -  <{»!.«?>  _  B1}/^2  -  <Bj,a2>  -1, 

L  <«  )  111  l^(Q  2 

p  p  P  p  Lp( 


8n  =  a"  -  <B^,ai?>Bi/ta”  -  \  2 

i=1  i=1  LJ 


n.a  ) 

P  P 


n-1 


For  this  scheme  to  be  well  defined,  we  need  to  show  that  1®^  -  £  ^  ( 

i=1  ' 


*  0.  Note  that  the  space  Sn  spanned  by  the  01,  i=1,***,n  is  the  same 
as  that  spanned  by  a1,  i  =  1,***,n,  because  the  matrix  An  =  [ai;j]  that 

>\  (e\\ 

carries  the  column  vector  j  J  J  to  I  J  I  is  lower  triangular. 


/*" 

I  •  •  • 


A  * 
n 


V 


\ 


J 


and  its  determinant  is  given  by  the  product  of  the  entries  in  the  diagonal,  a 


typical  entry  in  the  diagonal  is 

*u  -  *••}  -  rwj 


We  now  show  that 
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»«"  -  l  p  *  0,  Vn. 

i*  1  ' 


(6.  10) 


First  note  that  8*  =  — -  is  well  defined,  since  by  assumption  0  < 

la1 » 

1  2 ,  P 

la' I  <  •»,  Now  let  n  be  the  first  integer  greater  than  one  for  which  (6.10) 


fails,  and 


Then, 


<  -  l  • 

1*1 

-  X  sl 

n-1  . 

i  n  r\ .  r  « 1  n  .  1  a1  . 

’  (VV  ‘  ^  <P1'V2.P(W 
1=1 

-  («»  -  T  ^Xh.p  »;«"  -  "f 


-  (0.  a”  -  ”f  <eX>j,p  «2>- 
1*1 
n-1 

Since  an  e  W  and  I  „  6 e  W,  their  difference, 

i=1  ,P 

(o.  -  V  <»}.«?>,,»  »J>  «  ». 

However,  this  is  impossible  unless 

-2  -  j’  <6X>2,P  >2  *  °' 

For  if  5*0,  then  (5, A)  e  W  implies 

n  *  A  *  (-n  •  5>BV  on  T  , 

P 


n  *  A  *  0  on  T  . 

P 

But  as  mentioned  at  the  end  of  §4,  the  only  solution  on 

V  x  7  x  a  =  0, 

n  *  A  =  0  on  T  and  T  , 

P  v 


A  =  0. 
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Thus  we  have  shown  that 
n-1 


implies 


■?  - 1  * 

1=1 

“n  •  %  <•}•❖,., Bi  -  2  «6x>2.p  ?  -ji"1’ 

n-1  n-1  n-1 

‘  i,  'J,  <6X>2.p  V*  ■  j,  v  • 


n  n- 1  1 

In  other  words,  a  ,a  ,•••,<*  are  linearly  dependent  which  is  a 
contradiction. 

Now  since  as  mentioned  above  the  space  Sn  spanned  by  the  first  n 

vectors  an  is  the  same  as  that  spanned  by  the  0^,  the  closure  in  the  H 

norm  of  the  space  spanned  by  the  0  is  again  H.  Thus  assertion  (1)  is 

n 

proved. 

(2)  The  closure  of  the  set  of  all  finite  linear  combinations  of  0n  in  the 

W  norm  is  identical  to  that  of  o  and  hence  yields  the  space  W. 

n 

In  the  following,  we  define  the  Galerkin  approximations  to  the  solution 
of  the  AEVP.  Let 


Am(t)  =  (?m(t),Am{t))  =  l  c®(t)0i. 


We  shall  call  A  ( t)  a  Galerkin  approximation  of  the  order  to  the  AEVP 


-r  <  l  c"(t)P*,*J  >  +  a(<  l  c"*(t)0i),0j)  =  0 

t  i=1  1  1  1  ^/P  isl1  i 


V  0J  j  =  1, 


and  if  A  (t)  satisfies  the  initial  conditions, 

(0)  *  (£  (0),A  (0))  =*  (C0#aq) 

H,a.  tm 

H(0)  “  5o  ' 

Here  by  (6.5)  and  (6.6), 


(6.11) 
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<vv  ■  i  cX'62’ 

i=1 

00 

~0  „i 


-  a  3 »; 


i=1 


in  W, 


in  (H), 

P 


so  that 


<?“.»■>  -  i  f"  -  i  ct«i 

i=  1  i= 1 

Also  (6.11)  may  be  written  as 


dt 


/C1 


(t) 


+  B 


m1 


V 


=  o  , 


(6.13) 


m,  .  . 

C  (t) 

m 


where  B  =  [atS^/P^)]  is  an  m  *  m  matrix.  Note  that  the  coefficient  of 
m 

d2  i 

— —  C  (t)  is  the  identity  matrix  because  of  the  orthonormality  of  the  8 

d2t  "*  1 


in  L  (A  ). 
P 


The  theory  of  ordinary  differential  equations  now  ensures  a  unique 

solution  to  (6.13)  subject  to  the  initial  conditions  (6.12)  on  the  interval 

[0,T].  We  now  make  use  of  the  equations  (6.1)  to  obtain  a  priori  estimates  on 

l£m(t)l  and  IAm(t)l  .  Multiply  each  of  the  equations  (6.11)  by  cm(t) 

2 1 P  W  i 

and  summing  over  i  ,  we  obtain: 

<5m(t),£m(t)>2  P  +  a(Am(t) , Am(t) )  =  0.  (6.14) 

Now  because  of  the  symmetry  of  a( (C,A) , (5, A)  on  W,  the  relationship  may  be 
written 


4r  (iem(t)l2  „  +  a( Am(t) , Am(t) ) }  -  0.  (6.15) 

dt  <  fP 

Thus  after  integration  we  obtain 

l£m(t)l2  n  +  a<Am(t),Am(t))  -  lCm(0)l2  +  a ( Am ( 0 ) , Ara ( 0 ) ) .  (6.16) 

Since  the  bilinear  form  a( (C,A) , ( C,A) )  is  continuous  on  w,  there  exists 
c^  >  0  such  that 

|a(A“<0),Am(0))|  <  c1IAm(0)l2. 
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Also  by  virtue  of  Bessel's  inequality 

'^O'a.p  «  a<0),2,P  "  V 


Hence 


«Cm<t)«2  p  +  a(Am(t) , Am(t) )  <  c2  +  c^A^O)!2  . 


(6.17) 


Furthermore,  Am(0)  converges  to  A ( 0 )  in  the  W  norm,  so  there  exists 


c-j  >  0  such  that 
.m. 


Hence, 


I A  ( 0 ) I  <  c_.  Vm 
W  3 


,|m(t  ’  I  2  p  +  a ( Am { t ) ,  Am ( t ) )  <  c2  +  Clc3  -  br 


(6.18) 


Now  because  a(  (t,A) (C,A) )  is  coericive  with  respect  to  the  W  norm,  we  have 


a<  (£ ,A) , (€ .A) )  ♦  AKI2  p  >  6I(£,a)I2  , 


(t,A)  e  W. 

Thus  adding  vAl£l2  p  to  both  sides  of  (6.18)  and  using  the  coerciveness 


yield 


^m(t),2,P  *  6,<*m,Am>l£  «  ♦  Xl^m|l2,p 


It  follows  that 

Km(t)l2  <  b}  ♦  AKn,,2  p  . 

Multiplying  both  sides  of  this  inequality  by  I5ra(t)l2 

LP(£1P> 

Tt  «€"(t)f«-(t)>2  )  -  2«m<t),  ^  em(t)>2>p 


(6.19) 

(6.20 ) 

and  using 


*  2,tm(t),2,P1'^  5m(  ^  '  2  f  p  . 


we  obtain 


( V2^rK®<t)i2  J2  <  b.i^tt)!2  ♦  mm(t)i4. 

at  * ,  P  i  i  i  v 


Therefore, 


Since 


£-Um(t)l2.  <  2{b  ICra(t)l2  ♦  AKm(t,.4  ). 

dt  l2(Q  }  1  2,P  2,P 

P  p 
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we  obtain 


A»?r,m(t)ii*  +  b  ucm( t)  It^  <  < / X n cm ( t )  II J  + 

2  ,  P  I  ^  /  P  ^  j  P 


1  ,2 


2/A 


<  (/a  i^ra(t)»‘s  +  ~lr)2, 

/A 


”#4m(t)l!2  <  2(/A  ll^m(t)li2  +  — )  . 

at  2  ,  P  2  /  P 


Therefore, 


Km(t)« 


<  <>5*<0>l*  ♦ 

L„(Q  ) 

p  p 


b1  2 /At  b1 


Let 


b2  =  “^(0)"2,P  ' 
b2  =  supfb™}  . 

Note  that  b2  <  °°  is  ensured  from  the  inequality 

»5m(0)ll2  <  M4m(0)),Am(0)ll2  =  II  Am(0  )  II  2  . 

A  f  P  W  W 

Thus  given  e  >  0,  there  exists  M(e)  such  that  if  m  >  M(e), 
•An(0)l*  <  I A(0 ) «2  +  e  =  c,  +  e  , 

W  W  J 

because  Am(0)  +  A(0)  in  the  W  norm. 

Hence,  for  m  >  M(e) 

and  in  particular, 

=  lim  supt£m(0 )  II  2  <  c^. 

It  follows  from  (6.22)  that 

nl  y  1^4 

lim  suplC  ( t )  B 2  p  <  (c^  +  ~)e  -  j-  on  [0,T]. 


(6.21  ) 


(6.22 ) 


(6.23) 


Returning  now  to  equation  (6.19),  we  may  use  (6.23)  to  derive  similar  upper 


bounds  for  >^m(t)l2 


_  and  II  ( £m,Am)  (t )  II 2.  It  follows  that 

l2(«  )  W 

P  p 


lim  sup{  II  £m(t )  II  2  +  All  (£m,Am)  ( t )  U2  }  <  lim  sup{b  +  AIIE,mlt“  } 

2 ,  P  W  1  2 ,  P 


<  b  t  A(l5m(0)#2  +  -r1)e2,/H  -  b  =  A  (  #  ( 0  )  II 2  +  —  le2^.  (6.24) 

1  2 , PA  1  2  ,  p  A 


-46- 


From  (6.13)  and  (6.24)  together  with  e  <  e  on  [0,T],  we  conclude 


that  there  exists  M  >  0  such  that  on  [0,T]  ^ 

<  M,  Vm  , 

^  #  K 

^m(t)l2,P  <  Vm  ' 

M€m/Am)«f  <  M,  Vm  . 

W 

Thus, 

{ Ara ( t ) }  =  ((5ffl(t),Am(t)))} 
forms  a  bounded  sequence  in  the  Banach  space 

o>  2 

L  ( [0,T] :W)  cl  <{0,T]:W)? 

{£m(t)}  forms  a  bounded  sequence  on  the  Banach  space 

»  2 

L  ( [0  ,T] :W)  c  L  (  [0,T] ;H) . 

Therefore  we  may  extract  a  subsequence  such  that 
n  „ 

A  (t)  ♦  A(t)  in  L  ( [0,T] :W) , 

n 

A  m(t)  ♦  A(t)  in  LZ([0,T):W),  (6.26) 

•  n  ^  a, 

l  "(t)  ♦  ?(t)  in  L  ([0,T)]:Pfl  (H)). 

P 

In  particular  from  (6.26)  it  follows  that 

n  05 

C  m(t)  +  Ut)  in  L(  [0,T]  :LZ(S)  )). 

P  p 

Thus  £(t)  =  5(t),  because  of  the  uniqueness  of  the  representation  of  a 
vector-valued  distributional  derivative  (Lions  and  Magenes,  1972)  as  an 
element  of 

(l"([0,T1sL*(Q  ))  c  l\ [0,T]  :L*(0  )). 

P  p  P  P 

We  next  show  that  the  function  A(t)  defined  in  this  manner  is  indeed  a 
solution  of  the  AEVP.  For  this,  let  ^  e  CQ[0,T]  be  such  that 
i|»(0)  *  <J>(T)  *  0*  Let  'j'i(t)  «  <|> ( t ) B-*- ,  'l'j(t)  *  ♦( t ) We  then  multiply 
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(6*11)  for  A  by  ^(t)  and  then  integrate  by  parts  with  respect  to  t 


over  the  interval  [0,T]  to  obtain 
n  .  .  n 


/T< I  m(t),5^<t)>  +  a{ A  m(t),<l»1(t))dt  =  0. 

0  1  2,P 

Taking  the  limit  as  n  +  •  and  using  the  weak  convergence  of 

m 

n 

A  m(t)  to  A(t)  in  L2([0,T]:W),  we  obtain 

0  =  ft  <l(t),^(t)>2<p  +  a(A(t),K'1(t))dt 

*  JT{<l(t),BS  *(t)  +  a(A(t),B1)\)»(t)}dt 
o  1  2,p 

=  fT  t“T  4(t),B*>,  „  i(t)  +  a(A(t),B1)«|»(t)}dt. 

n  dt  l  2 ,  p 


(6.27) 


Thus 


/T  .  ♦(t)dt  =  J1  a(A(t)  ,Bi)'|»(t)dt 

n  dt  \  £,P  n 


0  —  ■  0 
for  all  *(t)  e  D[ 0/T)  . 
i  2 

We  show  g(t)  *  a(A(t),B  )  e  L  [0,T].  This  is  immediate  if  one  uses  the 
continuity  of  a(»,*)  in  the  W  norm,  i.e.  there  exist  c  >  0  such  that 
|a(A(t)  ,BA)  |  <  clA(t) I^IB1^  , 
and  if  we  pass  to  the  limit  in  (6.25),  we  obtain 
IA(t)«w  <  M. 


Therefore, 


so  that 


|a(A(t),Bi) |  <  CMlBilw,(6.28) 


g(t)  e  L  [0 ,T]  c  L  (0 ,T] 


(6.29) 


Hence  from  (6.28)  we  may  deduce 

a2  4  ? 

— r  <S(t),B,>,  e  L  (0,T) . 

dt2  *P 

The  derivative  is  again  taken  in  the  sense  of  distributions.  Furthermore, 

d2  i  i 

2  =  -a<A<t),e1)  . 

dt 
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form  a 


Now  since  this  equation  holds  for  arbitrary  01  and  since  the  01 

/v 

complete  set  in  W,  it  holds  for  every  (£,A)  e  W.  We  conclude  that 
d2 

~r  <£(t),5>.  +  a(A(t), (5, A) )  =  0  ,  (6.30) 

dt^  ' 

V(£, A)  e  w  , 

which  completes  the  demonstration  of  existence  of  a  solution  to  the  AEVP. 

§6.2.  Solution  of  the  EVP 

The  solution  to  the  EVP  is  constructed  in  the  following  way.  To  the 

solution  of  the  AEVP  defined  by  (6.5)  to  (6.9),  we  add  a  pair  (0,A*)  defined 

on  ft  *  ft  with  A*  being  a  solution  of 
p  v 

V  x  V  x  A* (t)  =  0  , 

V  •  A' ( t)  =  0  ,  (6.31) 

nXA'=0  on  T  u  F  , 

P  v 

/  (A'  *  n)ds  =  o(t)  on  [0,T], 

r 

v 

where  A'(t)  is  taken  from  V3  (see  4.4).  That  is, 

V  x  A'(t)  =0  n  x  A'(t)  =  0  on  T  u  T  , 

P  v 

V  •  A'(t)  =  0  /  A  •  n  ds  =  <J(t)  on  [0,T]. 

r 

v 

We  note,  in  particular,  that  if  we  define  A'  as  the  unique  solution  of 

V  x  A;  =  0, 

V  •  A'  =  0  , 

/  A'  •  n  ds  =  1  , 

r  1 

v 

A 

then,  A'(t)  *  <J(t)A^,  by  linearity  and  uniqueness.  The  solution  (£,A)(t) 
of  the  EVP  defined  by 

A 

(5,A)(t)  =  (S,A)(t)  +  (0,A')(t)  =  (5,A  +  A')(t) 
has  the  property 

A 

»(C,A)IW  =  »(C,A)«W  . 
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This  is  shown  as  follows.  Since  the  first  term  £  is  unchanged,  all  terms 
entering  into  1{£,A)HW  involving  only  £  are  unchanged.  The  vacuum 
contribution  is  given  by 

/  |7  x  (a  +  A')|2dv  =/  1 7  x  a  +  V  x  ft'|2dv 

ft  ft 

v  v 


since 


=  /  I  7  x  ft|2dv, 
ft 


7  x  A' ( t)  =  0 . 

Hence  we  have  (6.1) 

.2  .  _  ‘ 

— r  <£,£>_  n  +  a(  (£,A) ,  (£,A) )  =  0 
dt^ 

v(£,a)  e  w+  , 


and  also  the  natural  boundary  condition  L(£,A)  =  0  holds  whenever  the 
integration  by  parts  is  permissible  (§4).  The  initial  conditions  are  clearly 
also  satisfied  because  of  (6.6)  to  (6.8). 

In  summary,  the  solution  to  the  EVP  was  obtained  in  two  stages.  First  an 
AEVP  was  solved,  and  then  we  added  to  the  vacuum  part  of  this  solution  a 
solution  to  the  elliptic  system  (6.31).  As  noted  above,  this  elliptic  system 
with  the  given  flux  and  initial  conditions  has  a  unique  solution,  therefore  to 
show  the  uniqueness  of  the  solution  to  the  EVP  it  suffices  to  establish  the 
uniqueness  of  the  solution  to  the  AEVP. 

Let  A(t)  =  (£(t),A(t))  be  a  solution  to  the  AEVP  subject  to  the  initial 
conditions 


A(0 )  =  0  ,  £(0)  =  0. 

Let  s  e  (0,T).  Let  (£(t),A(t) )  be  defined  by 


£(t) 


-  /s  £ ( 0 ) do  t  < 
t 


t  >  s  , 
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A(t)  = 


-  J3  A(0)d a  t  < 


t  >  s  . 


Note  that  A(t)  =  (C(t),A(t))  e  W  Vt  6  [0,TJ,  simply  by  interchanging 
integration  with  respect  to  time  and  differentiation  in  space.  It  follows 


3  ^  3  ^ 

g^  £(t)  =  C(t),  g^  A(t)  =  A(t)  , 

A  A  A 

?(t)  =  Alt)  =  0  ,  t  >  s  . 
a2 

3  * 

Moreover,  note  that  — —  C ( t )  is  an  element  of  W  ,  the  dual  of  W.  Let 

3t 

32  ~  _ 

<_  2  *  denote  the  inner  product  between  W  and  W*  '*’n  t*'6 

3t 

duality  between  W  and  W*. 

Since  A(t)  is  a  solution  to  AEVP  with  zero  initial  conditions,  we  have 
•>2 


< — T  €,£>  *  +  a(A(t),A(t))  *  0. 

3t  w*w 


Hence, 


/t{<^-r  ?,£>  *  +  a( A(t) , A( t > ) }dt  =  0. 

0  3t  W*W 

♦ 

Integrating  by  parts  and  using  the  initial  condition  £{0)  =  0  yields 

JT{-<l(t),Ut)>  +  a(  A(t)  ,A(t) )  }dt  =  0. 

0  2,p 

Therefore,  since  (€(t),A(t))  =  (0,0)  for  t  >  s 

/S{-4(t),5(t)>  +  a(A(t),A(t))>dt  -  0, 

0  'v 

3  *v 

and  using  A(t)  *  A(t),  we  obtain 

/S  (|r  (K(t)l2  n  +  a(A(t),A(t))]}  =  0. 

0  3t  2fP 
It  follows  that 


a(A(0 ) ,A(0 ) )  +  IC(s)l  =  0. 


Now  adding  Al£(0)l  to  both  sides  and  using  the  coerciveness  of  a(*,»)  we 


obtain 


6IA(0)I^  +  K(s )*l  n  <  Al?(0)l2  , 

W  2 ,  P 
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or 


Thus, 


or 


IA(0) 1^  +  I5(s)lj  „  <  dl£(0)l2, 
W  2  ,P 


where  d  >  0.  We  set  W  =  (W^,W2)  with 


W  (t)  =  J  5(o)d<J  , 
0 

W  (t)  =  Jt  A ( a ) do . 

0 


«W(s)l2  +  *X(s)ll  <  dlw  ( s )  i2  , 

W  2.  f  P  1  2  /  P 


•w( S ) I2  +  l?(s)l2  p  <  dl/S  5(0)dal^  p  . 

rs  2 

Adding  the  positive  quantity  J  dlW(o)l  do  to  the  right  hand  side  of  this 

0  W 

inequality  we  obtain 


^  4>(s)  <  d*(s). 


(6.32) 


where  <f>(s)  =  /S  {lw(a)l2  +  l£(o)))l2  }dO.  Therefore,  since  <J> ( 0 )  =  0  , 
0  W  2,P 
integration  of  (6.32)  yields 

$(s)  =  0  and  -^  =  0  on  [0,T]. 
as 

Consequently , 


l£(t)l2  p  =  0  implies  5(t)  =0  on  [0,T]. 

It  follows  that  A(t)  =  0  by  Theorem  4.2.  Hence,  A(t)  =  (C(t),A(t))  =  0. 
This  establishes  that  the  only  solution  to  the  AEVP  is  the  null  solution. 

The  uniqueness  result  for  the  EVP  then  follows  trivially. 

In  conclusion,  we  note  that  the  solution  to  the  EVP  exists  for  all  time 
t  >  0.  This  result  follows  trivially  from  the  global  estimates  in  §6. 
Following  a  theorem  of  Lions  and  Magenes  (1972),  we  may  summarize  our  results 
as 

Theorem  6.1.  There  exists  a  unique  solution  to  the  EVP  (€,A)(t)  in 

C((0,*);w+)  with  £<t)  in  C((0,")jPjj  (H+))  such  that 

o 
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T  <5<t),|>  +  a( (C,A)(t),(£,A))  =0  V  (C,A)  e  w+, 

dt  ' P 


subject  to 


<S(0),A(0))  =  (50,aq)  e  w  , 


iLLQi  .  ?  e  p  (H+s 
at  S  e  pn  <H  '• 


and  the  flux  condition 


/r  A(t)  •  n  dS  =  o{t). 
v 


where  o(t)  is  a  prescribed  continuous  function  of  t. 
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§7.  Discussion 


First  we  make  some  remarks  about  the  solutions  to  the  EVP.  In  the 
process  of  establishing  existence  of  a  solution  to  the  AEVP,  we  construct 
Galerkin  approximation-  and  extracted  a  subsequence  of  the  approximations 
which  converges  to  the  solution  weakly.  We  note  that  the  entire  Galerkin 
sequence  is  in  fact  weakly  convergent  to  the  solution.  The  proof  is  based 
upon  a  well-known  compactness  property  (Yosida,  1978)  which  we  formulate  as 
the  following: 

Lemma; 

In  a  Banach  space  B,  the  first  two  assertions  imply  the  third  and 
conversely 


(1)  un  e  B,  un  weakly  compact, 

( 2 )  all  weakly  convergent  subsequences  of 
weakly  to  the  same  limit,  denoted  u. 


(3) 


un  is  weakly  convergent. 


We  make  use  of  this  lemma  in  the  following  way: 

oo 

Recall  that  A  was  weakly  compact  in  L  ( [0,T] :W) 
n 

•  oo 

(see  §6).  Also  5  was  weakly  compact  in  L  ( [0,T] :H) 
Consider  the  space 


un  must  converge 
imply  that 


2 

and  in  L  ( [0,T] : W) 
and  in  L2( [0,T] :H). 


B  =  L2(0,T]:W)  x  L2([0,T]:H). 

Let  un  *  (A^,^).  un  is  weakly  compact  in  B,  which  is  a  Hilbert  space 
with  norm 


MA  ,5  )l  =  {IA  «2  +  1 1  I2 

n  nL([0,T]:W)  n  L  ( [0,T] :W) 


Furthermore,  if  u  is  a  subsequence  which  converges  weakly  in  B,  then 

m 

exactly  the  same  sequence  of  arguments  as  given  for  (6.27)  to  (6.30)  shows 


that  the  limit  of  u  is  indeed  a  solution  to  the  auxiliary  evolutionary 
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variational  problem.  This  shows  that  the  entire  Galerkin  approximation  is 
weakly  convergent  to  the  solution. 

In  order  to  be  sure  that  our  solution  of  (EVP)  is  a  "good"  solution  of 
(MIBVP)  we  should  verify  that  it  satisfies  the  following  two  criteria: 

(1)  Any  classical  solution  of  MIBVP,  if  it  exists,  is  a  solution  of 
EVP. 


(2)  Any  classical  solution  of  EVP  is  a  classical  solution  of  MIBVP. 


To  verify  (1)  we  note  that  the  sequence  of  steps  leading  from  (3.1)  to 

<v 

(3.10)  performed  now  for  pairs  (C,A)  e  W+  shows  that 

•  •  A 

<P£,S>,  +  a( (€  ,A) , ( C,A) )  =  /  (?  •  n)L( £,A) ds, 

^ ,  P  P 

P 

if  (€,A)(t)  is  a  classical  solution  of  MIBVP.  The  right  hand  side  however 


vanishes,  since  the  boundary  condition  (3.4)  is  satisfied  by  a  solution  to 


MIBVP,  if  it  exists. 


To  verify  (2)  we  start  from 

g2 

— -  <S,£>  „  +  a((C,A),(?,A))  =  0, 

3t2  2,P 

^  ~ 

(£,a)  e  w+, 

and  assume  (£,A)(t)  is  a  C2  function  on  ft  x  ft  ,  and  C2  in  time. 

P  v 

We  then  may  integrate  by  parts  the  various  integrals  in  a( (S,A) , (5, A) ) , 
d2 

and  bring  the  operator  — r  inside  the  scalar  product  to  get 

dt2 

2  •v 

—  <5(t),C>_  „  +  a((C,A),(C,A))  = 

A*.* 


"  C(t),5>2>p  -  <F(C),?>2fP  t  <V  x  V  x  A,A>2  v 


=■  /  (5  •  n)L(C,A)ds, 

r 

P 

where  Aew+,  VxVxa=o. 
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W  ^  w 

furthermore  since  C„(ft  )  Pn  (W  ),  we  may  choose  for  4,  any  Cn(ft  ) 

up**  up 


function,  so  that  the  boundary  contribution  vanishes,  therefore 

32  2 

— r  4(t)  =  F(4(t))  in  D*([0,T] :C  (ft  )>, 

8t2  p 

(Brezis,  1974),  since  F(4(t))  e  C°[0,T]:C2(ft  )], 

P 

4(t)  e  C2U0,T]  :C2(ft  )). 

P 

13+  ^ 

Finally,  since  H  (ft  )  Pjj  (W  ) ,  we  may  take  for  (4  •  n)  any  element 

2  P  P 

in  L  (T  )  so  that 
P 

L(4,A)  =  0  in  L2(T  )  =  (L2<r  ))*. 

P  P 

Thus,  since  1,(4, A)  is  represented  by  an  element  in  C^(T  ), 

P 

1.(4, A)  =0  in  C°(T  ). 

P 
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APPENDIX 


We  now  establish  the  identity 

7(5  •  7  x  b  *  B)  =  7(5  •  7b)  x  b  +  7  x  b  x  (5.  7b) 

-  [ ( (75 )T  •  7)  x  B]  x  b  +  (7  x  B  x  B)  •  75.  (A. 1 ) 

This  identity  need  not  be  true  for  all  pairs  5  and  B,  and  only  for 
those  where  B  satisfies: 

7  •  B  =  0 

V  x  b  x  b  =  7p  (A. 2) 


and  we  assume  that 

(7  x  b  x  b)1  =  (7  x  b  x  b)3  (A. 3) 

J  1 

which  would  follow  if  p  is  twice  continuously  differentiable  (equality  of 
the  mixed  second  partial  derivatives).  To  prove  the  identity  we  first 
establish  its  validity  for  5=6^  5  =  or  5  *  e^  the  Cartesian  basis 

vectors  and  B  satisfying  (A. 1). 

For  5  *  e^,  the  identity 

7(5  »7xbxb)=7x(5»  7b)  xb+7xbx(5»  7b) 

-  ((75)t  »7)xbxb+7xb*B, 


becomes 


since 


7(7  x  b  x  b)  *  7  x  ( 


3x, 


B)xb+7xbx( 


3x. 


B), 


1 

Now  the  R.  H.  S.  can  be  written 
3  (7  x  b  x  b). 


so  the  identity  (A.1)  reduces  to  (A. 3). 

51  and  52  it  is  clearly  true  for  5 
a  sum  of  terms  involving  51  and  5?» 
it  is  true  for  5  it  is  true  for  o(x) 


Now  if  the  identity  is  true  for 
*  51  +  52  since  all  terms  split  into 
it  is  slightly  harder  to  show  that  if 
where  a(x)  is  a  function  of 
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i.e.  we  seek  to  show: 


x, ,x2,x3, 

V(a£  *V*BXB)  =  7  x  ((a£)  •  7b)  *  B 

+  V  x  B  x  ((aC)  •  7b)  -  (7< aC)T  •  7)  x  B  x 
+  V  x  b  x  b  •  7(a£), 

follows  from 

V(?  •  7  x  B  x  B)  =  7  x  (C  •  7B)  x  B  +  7  x  B  x  (5 

-  <(7?)T  *  V)xbxb+7xbxB*  75. 


For  this  note  that  using 


7 ( fg )  *=  (7f)g  +  <7g)f, 

the  L.  H.  S.  of  (A. 4)  is 

a7(5  •  7  x  b  x  b)  +  (7a)((7  x  b  x  b)  *  5). 


Using  (A. 6)  and  the  identity 

7  x  (fv)  =  f7  x  v  +  7f  x  v  , 

the  R.  H.  S.  becomes 


i(7  x  (5  •  7b)  x  b)  +  {7a  x  (5  •  7b)}  *  B 

+  a( (7  x  b)  x  (5  •  7)B  x  5  -  a[ (75)T  •  7]  x 


-  [  |ayj  ’  7J  *  B  *  B  +  *  B  * 

+  (7  x  b  x  b  •  £)7a. 


where 


a 

x 


3a 
3x  ' 


a  «  -5— 

y  3y 

a  »  -s— 
z  3z 


Cancelling  (7a)  (7  x  B  x  b  •  5)  from  both  sides  and  using  (A 


to  show 


{7a  x  (5  •  7b)}  x  b 


U 


/ax\ 
I  ay 
\az  J 


7} 


B 

(A. 4) 

♦  7B) 

(A. 5) 

(A. 6) 

B  x  b 

B  •  75 

.5),  it  remains 

B  x  b. 
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In  fact. 


The  L 


and  the  R. 


which  are 


1 


ax  ’ 


7a  *  (£  •  Vb)  =  t  [  ay)  *  71  *  B* 

H.  S.  is 


,  az , 


det 


3  _i 


I  5lieB 


ay 

iT3  C1  3  B2 
i*1  i 


k 

az 

i=3  ■  a 

l 

i=l  i 


H.  S.  is: 


detj  1  ax  tX  g^- 
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